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A simultaneous description of non-strange nuclei, hypernuclei and multiply-strange nuclear sys- 
tems is provided by a single mass formula which is shown to be useful for estimating binding energies 
of nuclear systems over a wide mass range, including the light mass nuclei. It not only provides 
a good fit to the existing experimental data on hyperon-separation energies but also reproduces 
■ results of the relativistic mean field (RMF) calculations. In addition, it can provide the Lambda(A), 

' Cascade-0(H°) and Cascade-minus (S~) drip lines. The existence of a range of bound pure-hyperonic 

systems without any neutron and proton is suggested among which 6A, 9H°H°, lOSTST, 1A7H , 
1A8E - , 1H°9H~, 1H~8H° and 2A + 3H + 3H~ represent the lightest species. In agreement with the 
RMF predictions, this generalized mass formula also predicts the nucleus 2s°2A 8 He to be bound. 
An exotic 2H°2H-2A 10n nucleus is also found to be bound. This new mass formula can be used in as- 
, trophysics for strange stellar objects, as well as in high energy physics for estimating the strangeness 

^s^j ■ production yield in nucleus-nucleus or, nucleon-nucleon collisions. 
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Hyperons such as Lambda(A), Cascade (2) are unstable and they lose their strangeness by decaying via the weak 
interaction to lighter particles containing only up and down quarks. Bodmer [l[ and Witten Q hypothesized that 
for a large collection of quarks a roughly equal number of up, down and strange quarks leading to a strangeness 
fraction f s — S/A w 1 and charge fraction f q — Z/A k, would form a stable state of strange quark matter 
|^ (SQM) called, "strangelet" @. Strangelets are conjectured as a possible candidate for clean energy source, and in 
astrophysics, for dark matter. Strangelet search at high energy heavy ion collision e.g, 200 GeV Au+Au Collisions 
at the BNL relativistic heavy ion collider set upper limits for strangelets with mass > 30 GeV/c 2 Strangelet pro- 
duction possibility at the Large Hadron Collider (LHC) has been investigated @ and strangelets are yet to be detected. 

rn ; 

In nineties, Schaffncr ct al. pointed out that at 2-3 times the nuclear matter density, metastable strange hadronic 
matter (SHM) made of Lambda (A) and Cascade (S) hyperons along with nucleons (N) might also exist with 
f s = S/A ss 1 and f q = Z/A « @,0]. Like strangelets, metastable multi-strange hypernuclei could also be produced 
in heavy ion collisions 8, 9j . The metastability of the strange hadronic systems was established @, 0] by extending 
relativistic mean field (RMF) calculations from ordinary nuclei (f s = 0) to multi-strange nuclei with (f s ^ 0). From 
a relativistic meson-baryon field theory metastable multi-hypernuclear objects and negatively charged composite 
objects with positive baryon numbers were predicted to exist [9]. In fact the RMF calculations have shown that if 
the hyperon- hyperon interaction is attractive, truly exotic pure-hyperonic nuclei (MEMO's) can exist d, S G3 • In 
1977, R. Jaffe had proposed existence of a possible stable AA called, "H dibaryon" (uuddss), made by combining two 
A (uds) hyperons [ll| and in 2006 G. A. Miller suggested existence of a possible double cascade (SS) di-baryon [l2l |. 
It has also been suggested that there could be bound hypernuclei whose counterparts in normal nuclei are unstable 
(such as, the experimentally found bound j^Li hypernucleus [13j | whose normal counterpart 10 Li is unbound) or, 
may even be beyond the normal neutron-dripline [LH [l5j . As the knowledge of hyperon- hypeon interaction or, 
hyperon-nucleon interaction is still not well understood, it is necessary to produce and study bound hypernuclei of 
different strangeness and mass numbers. 

In a mean field approach, obtained from a set of phenomenological meson-baryon interactions, Schaffncr ct al 
[3] did an extensive study on the stability of multiply-strange baryonic systems. Such systems were found to be 
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weakly bound (Eb/A «-10 to -20 MeV), but can be stable for large A. The stability mainly arises by the presence 
of S~ hyperons which reduce the Coulomb repulsion contribution in binding energy from protons. When a large 
amount of strangeness is incorporated into the system, the mean fields felt by a nucleon (N) and by a hyperon 
(Y) are altered, since YY as well as YN interactions become important. Schaffher et al [7] presented two models 
in the RMF framework; one without (Modell) and one with (Model2) the YY interactions. However, such RMF 
calculations are difficult to carry out for all kinds of nuclei, especially if they are isospin asymmetric. Immediately, 
Balberg et al [16] recognised that the use of a properly constructed mass formula would provide a quick check 
on the more complex RMF calculations and it would also allow to reliably extrapolate into domains not explored 
by RMF or, other methods. In 1993 Dover and Gal had proposed (starting from the basic (p, n) Fermi-gas 
model and then extending it to (p, n, A) and then (p, n, A, H°, 3~) systems) a generalised Bethe-Weizsacker mass 
formula which predicted a wide range of stable multi-strange nuclei. Balberg et al [l6] updated the parameters of 
the said mass formula in order to compare as close as possible with the results of the RMF calculations 0- They 
tried to reproduce the above RMF results by using two sets of parameters, Set-I and setll (see Table 1, Ref. 
corresponding to the RMF calculation's Modell and Model2 respectively. As this mass formula was used to reproduce 
the RMF results of multiply-strange nuclei, it was not tested to see how it reproduces the available experimental data 
on hypernuclei. Nevertheless, it demonstrated the usefulness of a generalised mass formula for strange hadronic matter. 

Relativistic heavy ion collisions provide a prolific source of strangeness. With the advent of high energy accelerators, 
the production of strange matter as well as dibaryons with strangeness have become topics of great current interest 
(l8j . Searches on a variety of strange-nuclei are currently being pursued at different laboratories around the world 
such as, JLab, J- PARC, GSI etc. The up-coming experiments at LHC (CERN) are also expected to produce 

some exotic strange matter. In this context, Botvina and Pochodzalla [2(| showed further usefulness of a mass formula 
of hypernuclei. They used the generalised mass formula for non-strange normal nuclei and strange hypernuclei 
of Samanta et al [lj], and the liquid drop model of W. Greiner [2l[ to calculate several quantities such as, the 
strangeness chemical potential versus temperature, average number of hyperons in fragments, and yield of fragments 
with more than one hyperon produced in a high energy heavy ion collision. Needless to say those calculations are 
extremely useful as knowledge of these quantities can provide guidance to the upcoming experiments looking for 
unknown strangeness production. 

In 2002, Samanta and Adhikari extended the Bethe-Weizsacker mass formula for light nuclei in which the binding 
energy of a nucleus of proton number z c , neutron number n and total baryon number A = z c + n is given as [22[ 

B(A, z c ) = 15.777A - 18.34A 2 / 3 - 0.71z c (z c - I) /A 1 ' 3 - 23.21(n - z c f /[(l + e~ A/17 )A] + (1 - e- A ' 30 )8 

(1) 

where 8 — — 12 A^ 1 / 2 for odd n-odd z c , — 12 A^ 1 / 2 for even n-even z c , and = for odd A. 

Later on, the above mass formula was extended to achieve a simultaneous description of non-strange nuclei and 
singly-strange hypernuclei in which the hypernucleus was considered as a core of normal nucleus plus the hyperon(s), 
and the binding energy was defined as [14j 

B(A, Z) = 15. 777,4 - 18.34A 2 / 3 - 0.71Z(Z - I)/ A 1 ' 3 - 23.21(n - z c ) 2 /[(l + e~ A / 17 )A] + (1 - e - A / 30 )5 
+ny[0.0335(m r ) - 26.7 - 48.7 | S \ /A 2/ % 

(2) 

where ny — number of hyperons of a particular kind in a nucleus, my = mass of that hyperon in MeV , S = strangeness 
of that particular hyperon and mass number A = n + z c + ny is equal to the total number of baryons. Note that the 
mass (m A = 1115.683 MeV, m B - = 1321.71 MeV, and to h o = 1314.86 MeV) and strangeness (S A = -1, S s o = -2 
and SW- = —2) of hyperons are explicitly considered in this formula. The total charge of the strange system can be 
negative if no- is more than z c . Therefore the Z in eqn.(2) is given by 



Z = \z c + n Y q Y \ (3) 

where qy is the charge number (with proper sign) of hyperon(s) constituting the hypernucleus. For non-strange 
(S=0) normal nuclei, z c — Z as ny =0. The choice of 5 value still depends on the number of neutrons and protons in 
both normal and hypernuclei. For example, in case of \Li the neutron number n=5(odd), proton number z c =3(odd), 
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and riy= 1. Therefore, 6 = -12^- 1 / 2 as the (n, z c ) combination is odd-odd, although the total baryon number 
A = n + z c + ny =9(odd). Whereas, for non-strange normal 9 Li nucleus 8 — for A=9(odd). 

The hyperon separation energy Sy was defined as [l4| 



Sy = B(A, Z) hyper - B(A - n Yl z c ) core (4) 

This mass formula produced good fit to both the normal nuclear binding energies as well as the available experimental 
data on binding energies of A, AA and S~ hypernuclei [14| . But, the binding energy per nucleon versus A plots 
delineate a significant discrepancy between the results of this mass formula and the RMF predictions [7j for nuclear 
systems with large A values. 

In this work the above mass formula for non-strange and singly-strange hypernuclei (i.e., with only one kind of 
hyperons)is extended for multiply-strange nuclear system which can have a mixture of Y= A, S° and H _ . The 
binding energy expression, found by a simultaneous reproduction of experimental data of binding energies of normal 
and strange hypernuclei along with the theoretical RMF calculations for multiply-strange nuclear systems, is given as 



B(A, Z) = 15.777,4 - I8.34.4 2/3 - 0.71 Z{Z - 1)/A 1/3 - 23.21(n - z c ) 2 /[(l + e - A/17 )A] + (1 - e _A/30 )<S 
+ n Y [0.0335(m Y ) - 26.7 - 48.7 | S \ /A 2/3 

Y 

-a Y {{n A + n E a + n s - - z c f + (n A + n S o + n s - - n) 2 }/{{\ + e - A/17 )A}}. 

(5) 

Here only one new parameter ay is introduced to consider hyperon-number asymmetry with respect to the proton 
and neutron numbers. This hyperon-proton and hyperon-neutron asymmetry term, introduced at the end, has 
been incorporated in a way similar to the neutron-proton asymmetry term of the normal non-strange nuclei 
(eqn.l). For ny=0, the hyperon terms disappear and the equation becomes that of normal nuclei. Addition of 
the hypeon-asymmetry term alleviates the above cited discrepancy at large A (Fig. la), while it does not have any 
significant effect on the results for lighter mass nuclei (Fig. la and Fig. lb). Interestingly a choice of the value ciy=0.02 
reproduces both the experimental data (Fig. 2) as well as the the RMF (Model 2) results over a wide mass range 
(Fig. la and Fig. 3) [tJ. Further RMF (Model 2) calculations for nuclei with large neutron-proton asymmetry would 
be useful for checking the choice of this ay value. 

The total charge of the multiply-strange nuclear system can be both positive or, negative for positive baryon 
numbers and in the above equation the charge Z is taken as 



Z = 



z c + 



/ nyqy 

Y 



(6) 



Note that q\ = g~o = 0, but q^- = —1- Therefore, the total charge of a hypernucleus depends on the number of 
protons (z c ) and the number of S - hyperons (n^- ), i.e., Z — \z c — n-^-\. This makes hyperonic systems with 
only S° hyperons more bound than pure S~ systems (as the negatively charged hyperons will encounter repulsive 
Coulomb force). On the contrary, for nuclear systems with both protons and hyperons, the coulomb repulsion of the 
negatively charged S~ hyperons are reduced by the protons, making S~ more bound than 3°, until the hyperon 
number exceeds the proton number. 

The Model 1 of the RMF calculations of Schaffner et al 7} can also be reproduced by adding a correction term C r 
to the above expression of binding energy. But, the origin of this term is not clear. Nevertheless, it is given in the 
following for the sake of completeness. We define 



B(A,Z)cs2 = B(A,Z) 
B{A,Z) CS i = B{A,Z)~C r 

(7) 
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The correction term C r and the strangeness fraction f s are given by 

C r = 12.0Af s (f A + f B o+f s -) 

(8) 

f s = Y,n Y \S\/A 

Y 

(9) 

The /a = n\/A, = n-~°/A, / s - = n s -/A and the sum of the proton fraction (f. p — z c /A), neutron fraction 
(fn = n/A) and hyperon fraction (/y = ny/A) makes /„ + f p + J2y f Y = 1- ^ n ^ nc following the results of this 
new mass formula are compared with the experimental data on hyperon separation energies, RMF calculations Q 
on binding energies of multiply-strange nuclear systems and calculations of Balberg et al In the figures the 

binding energy, B(A, Z)cs2, which reproduces the RMF(Model2) calculations with YY interaction, is represented as 
'CS2', and the other one (B(A, Z)csi) which reproduces RMF(Model 1) calculations without YY interaction [7| as 
'CS1'. The notation Eb/A = —B(A,Z)/A is used to compare the results of this new mass formula with the RMF 
calculations of Schaffher et al Q. 

Fig. 2 shows that both CS1 and CS2 give good fit to the experimental data. A comparison of the RMF calculations, 
Balberg's model2(Bal2) and modell (Ball) results with those of CS2 and CS1 are shown in Fig. 3 and Fig. 4. The 
most striking aspect of these plots is that how well the RMF predictions agree with the CS1 and CS2, especially in 
the low A value region. 

The Eb I A versus A plot (Fig. 5a) of pure hyperonic matter made of equal number of A, 5° and 5~ shows that at 
larger A the predictions of CS2 are different from that obtained using the formula of Balberg et al (Bal2), and the 
zero separation energy of A = 5° = 5" = 1 (Fig. 5b) are entirely different for Bal2 and CS2. The lowest bound pure 
hyperonic matter made of A, 5° and 5" (without any neutron and proton) is found to be 2A + 35° + 35~ (Eb /A 
= -0.603 McV) which is close to the value 2A + 25° + 25" predicted by the RMF calculations. The symmetric 
combination 3A + 35° + 35" (Eb/A = -2.947) is also bound. The upper limit for the A (=n A + n-=o + n H o) value 
of the of pure hyperonic matter (Fig. 5b) is (A ~ 285). In this mass formula, the H dibaryon is not bound. A 
bound system can be found with minimum 6A and maximum 304A hyperons. For Cascade, the 95° (Eb/A = -1.57 
MeV) and 7485° (Eb/A = -0.002 MeV) are the minimum and maximum numbers required for binding while for 
the 3~ the first bound is 105" (E B /A = -0.63 MeV) and the maximum is only 1555" (E B /A = -0.01 MeV). The 
1A + 75° (Eb/A =-0.10 MeV) is the minimum number of 5° that will bind with A but, 7A + 15° is more strongly 
bound (E b /A = -4.22 MeV). The 1A + 85" (E B /A = -0.16 MeV) is the lowest bound for A, 5" combination. The 
IS" + 85° (E B /A = -1.60) is bound and, 15" + 75° (E B /A = 0.563) is not. The other combination 85" + 15° 
is not bound (E B /A = 0.35 MeV) and the minimum number of S~ that will bind is 15° + 95" (E B /A = -1.20 
MeV). The nuclei 2 E<>2A 8 He (Eb /A = -4.51) and 2 e-2A 8 ™ ( A He + 25" + 2A, E B /A = -4.66) are also bound in 
agreement with the RMF predictions. The exotic 2~°2H-2A 10rl is found to be even more strongly bound(£ , s /A = -5.98). 

The usefulness of this mass formula is not only for the binding energy of normal as well as SHM, it can also give 
the hyperon drip lines. The generalised hyperon separation energy Sy is defined as 

Sy = B(A,Z) - B(A-^2n Y ,Z) (10) 

Y 

Note that the Z value for each nucleus has to be calculated considering its proton and 5~ hyperon numbers, and for 
nuclei wih A and 5° the total charge Z = z c . The fig. 6 and Fig. 7 show the hyperon-drip points (Sy = 0) for 56 Ni 
and 310 G cores plus hyperons predicted by both CS1 and CS2. The CS2, which reproduces the RMF results with YY 
interaction, yields a higher limit for hyperon separation energy (as expected) than CS1. The figure clearly indicates 
that there is a maximum value for the hyperon number which can be bound in a nucleus and that end point would 
differ according to the neutron and proton contents of the core nuclei. The higher the proton number, the larger 
the value of the drip point. Joining these drip points one can get the hyperon driplines for the entire nuclear chart. 
Fig. 8 represents the variation of binding energy per nucleon (Eb/A) with the hyperon fraction (f s ) for multi-strange 
hypernuclei with 56 Ni and 1S0 Th cores. In Fig. 8 a)-b) comparisons have been shown with a choice of A, 5° and S~, 
made in RMF calculations @,[H, CS2 and Bal2(with same combination of A, 5° and 5" as RMF), CS2(1,1,1) and 
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Bal2(l,l,l) in which the A, 5° and 5~ values have been increased by 1A, IS and IS - respectively. At larger f s 
here is a slight difference between CS2 and CS2(1,1,1); as well as Bal2 and Bal2(l,l,l) for the same core and it arises 
due to the choice of hyperon combinations. Both plots indicate an upper limit for hyperon fraction in finite nuclei. 

In summary, in this work we extend the hypernuclear mass formula of Samanta et al [14| for multiply-strange nuclear 
systems and present a generalised mass formula for non-strange, singly-strange and multiply-strange nuclear systems. 
It reproduces the RMF results of Schafiher et al [7| as well as the existing experimental data on hyperon-separation 
energies from hypernuclei. It suggests existence of a range of bound pure hyperonic matter of various combinations. 
The lighest combinations are 6A, 95°S°, lOS^S", 1A75 , 1A85", 1S°9S~, 1S"8S° and 2A + 35° + 33". The 
calculation shows that there are upper-limits for binding of the pure hyperonic matter as well as the strange and 
multi-strange nuclei leading to respective drip lines. The hyperon-drip points of hypernuclei with 56 Ni : 20S Pb and 310 G 
cores and A, 3 , S~ hyperons are presented. Because of its simplicity, this new generalised mass formula can be easily 
used in astrophysics (for strange stellar objects), as well as in relativistic heavy ion collisions for calculating different 
parameters such as, multiply-strange matter production yield, strangeness chemical potential versus temperature etc. 
as carried out [20} with the singly-strange mass formula of Samanta et al [lH and liquid drop model of Greiner 



21]. Such information is urgently needed for upcoming heavy-ion or, nucleon-nucleon collision experiments which are 



expected to produce multiply-strange nuclear systems. 
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FIG. 1: Effect of the new hyperon-asymmetry term (ay):(a) Eb/A versus A plots for stable multi-strange systems in RMF 
calculations based on 56 Ni nuclear cores for Model2 (with Y-Y interaction) of Schaffner et al. Q, and eqn.5 with ay = 0.0, 
0.015, 0.02 and 0.025 (b) Single Lambda-hyperon separation energy Sy versus mass number A plots with ay =0.0 and ay 
=0.02 for different elements. 
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FIG. 2: Hyperon separation energy Sy versus mass number A plots for (a) experimental values (exp) (tabulated in Ref. 
for single A (my = 1115.683 MeV) and double-Lambda (inset) predicted by Set-I (Ball) and Set-II (Bal2) of Balberg et al 
and CS1 and CS2 (This work); (b) the same for single S~(m s - = 1321.71 MeV) -hyperon. 
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FIG. 3: Eb/A versus A plots for stable multi-strange systems in RMF calculations based on 5e Ni and 180 Th nuclear cores for 
Model2 of Schaffner et al 0, Set-II of Balberg et al (Bal2) and this work (CS2). 
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FIG. 4: E B /A versus A plots for stable multi-strange systems in RMF calculations based on 56 iVi, 132 Sn, 208 Pb and 310 G 
(z c = 126, n = 184) nuclear cores for Modell of Schaffner et al 0, Set-I of Balberg et al (Ball) and this work (CS1). 
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FIG. 5: Plots for pure hyperonic multiply-strange systems consisting of baryon number A = ha + n s a +n 2 - (a) binding energy 
per nucleon Eb/A versus A and (b) separation energy Sy of 1A + IS + 1H~ hyperons versus A. 
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FIG. 6: Plots (a)-(c) are for A, 3°, and S separation energies versus A for A, H°, and 3 hypernuclei with 56 Ni core to 
estimate the drip point (5V = 0) and (d)-(f)are for binding energy versus A calculated by CS1 and CS2. 
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FIG. 7: Plot of A, S , andS separation energies versus A for A, S , and 5 hypernuclei with G and Pb core to estimate 
the drip point (5V = 0) of the respective nucleus. 
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FIG. 8: Eb/A versus f s plots for multi-strange hypernuclei for RMF, CS2 and Bal2 (with same hyperon combinations as RMF 
BUS]); CS2(1,1,1) and Bal2(l,l,l) in which the A, H° and E~ values have been increased by 1A, 1H° and IS - respectively; 
(a) with 56 Ni core (b) with 180 Th core. 



